Radiative transfer theory for vacuum fluctuations 
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A semiclassical kinetic theory is presented for the fluctuating plioton flux emitted by a disordered 
medium in thermal equilibrium. The kinetic equation is the optical analog of the Boltzmann- 
Langevin equation for electrons. Vacuum fluctuations of the electromagnetic field provide a new 
source of fiuctuations in the photon flux, over and above the fluctuations due to scattering. The 
kinetic theory in the diffusion approximation is applied to the super-Poissonian noise due to photon 
bunching and to the excess noise due to beating of incident radiation with the vacuum fluctuations. 
PACS: 42.50.Ar, 05.40.-a, 42.68.Ay, 78.45. +h 



The theory of radiative transfer was developed by 
Chandrasekhar ^ and Sobolev |^ to describe the scat- 
tering and absorption of electromagnetic radiation by in- 
terstellar matter. It has become widely used in the study 
of wave propagation in random media, with applications 
in medical imaging and seismic exploration |^ . The basic 
equation of radiative transfer theory is a kinetic equation 
of the Boltzmann type, that is derived from the Maxwell 
equations by neglecting interference effects It is a 
reliable approximation whenever the scattering and ab- 
sorption lengths are large compared to the wavelength, 
which applies to all but the most strongly disordered me- 
dia. 

Radiative transfer theory has so far been restricted 
to classical waves, excluding purely quantum mechani- 
cal effects of vacuum fluctuations. This limitation is felt 
strongly in connection with the recent activity on ran- 
dom lasers These are amplifying systems in which 
the feedback is provided by multiple scattering from dis- 
order rather than by mirrors, so that radiative transfer 
theory is an appropriate level of description. However, 
while stimulated emission has been incorporated into this 
approach a long time ago by Letokhov j^, spontaneous 
emission has not. It is the purpose of our work to remove 
this limitation, by presenting an extension of the radia- 
tive transfer equation that includes vacuum fluctuations 
and the associated spontaneous emission of radiation. 

Our inspiration came from the field of electronic con- 
duction in disordered metals, where the notion of a 
fluctuating Boltzmann equation (or Boltzmann-Langevin 
equation) has been developed extensively , following 
the original proposal by Kadomtsev . In that context 
the fluctuations originate from random scattering and 
they conserve the particle number. This same class of 
fluctuations exists also in the optical context considered 
here, but with a different correlator because of the differ- 
ence between boson and fermion statistics. In addition, 
the photons have a new class of fluctuations, without 
particle conservation, originating from random absorp- 
tion and emission events. Vacuum fluctuations are of 
the second class. We will extend the radiative transfer 
theory to include both classes of fluctuations. To demon- 
strate the validity of our "Boltzmann-Langevin equation 



for photons" , we solve the problem of the excess noise 
from vacuum fluctuations in a waveguide geometry, for 
which an independent solution is known We then 

apply it to the unsolved problem of the thermal radiation 
from a spherical random medium. 

The basic quantity of the kinetic theory is the fluc- 
tuating distribution function /k(r, i) of the number of 
photons per unit cell {2tt)^^ dkdr in phase space. (For 
simplicity, we ignore the polarization dependence.) Con- 
ventional radiative transfer theory deals with the mean 
/k(r), which we assume to be time-independent. It sat- 
isfles the Boltzmann equation 



■ ^ = E('^kk'(/) - Jk'k(/))+4+(/) - /k (/). (1) 



[For ease of notation, we write instead of (27r)"^ J dk, 
and (5kq instead of (27r)'^(5(k — q).] The left-hand-side 
is the convection term (with c the velocity of light in 
the medium and k a unit vector in the direction of the 
wavevector k). The right-hand-side contains gain and 
loss terms due to scattering, Jkk'(/) — fi'kk'/k'(l + /k), 



due to ampliflcation, I^{f) 



(1 + /k), and due 



to absorption, I^{f) = Wk/k- The scattering rate 
uikk' = Wk'k is elastic and symmetric. The absorption 
and ampliflcation rates are isotropic (dependent only 
on fc = |k|) and related to each other by the requirement 
that the Bose-Einstein function 



f,q{u,T)^[expihLu/kBT)-l]- 



(2) 



is the equilibrium solution of Eq. (|i|) (at frequency uj — ck 
and temperature T). This requirement fixes the ratio 
Wj^/wj^ = exp (TiLLi/kBT). The temperature T is positive 
for an absorbing medium and negative for an amplifying 
medium such as a laser p^ ] . 

We now extend the radiative transfer equation (^ to 
include the fluctuations 5f = f — f. Following the line 
of argument that leads to the Boltzmann-Langevin equa- 
tion for electrons ||7|-|l0[| , we propose the kinetic equation 



^k-^ = E(-^kk'(/)-Jk'k(,f)) 



(3) 
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The argument is that the fluctuating / is propagated, 
scattered, absorbed, and amphfied in the same way as 
the mean /, hence the same convection term and the 
same kernels Jkk',^k^ appear in Eqs. (|l|) and (|^). In ad- 
dition, Eq. (||) contains a stochastic source of photons, 



E 

k' 



kk' 



5J] 



k'k) 



k ' 



(4) 



consisting of separate contributions from scattering, am- 
plification, and absorption. This Langevin term has zero 
mean, — 0, and a correlator that follows from the as- 
sumption that the elementary stochastic processes (5 Jkk' , 



SI^ have independent Poisson distributions: 



5Jkk'(r,t)(5Jqq'(r',i') = A(5kq(5k'q' Jkk' (/), 

SI^{r,t)SI^{r\t') = AJkq/^ (/), 



5Jkk'(r,t)<5i^(r',t') = 0, «+(r,i)(5/q (r',t') = 0, 

where we have abbreviated A — 5{r — r')S{t — t'). 
stitution into Eq. (Eh gives the correlator 



(5a) 
(5b) 
(5c) 

Sub- 



£k(r,t)/:q(r',i') 

- ^kq(/) - 



<5kq^(Jkk'(/) + Jk'k(/)) 



Jqk(/)+<5kq(4+(/)+4"(/)) • (6) 



Eqs. (g) and (g) constitute the Boltzmann-Langevin 
equation for photons. 

To gain more insight into this kinetic equation we make 
the diffusion approximation, valid if the mean free path 
is the shortest length scale in the system (but still large 
compared to the wavelength). The diffusion approxi- 
mation consists in an expansion with respect to k in 
spherical harmonics, keeping only the first two terms; 
/k = /o + k • fi, /Ik = -Co + k • £1, where /o, fi, £0, and 
Ci do not depend on the direction k of the wavevec- 
tor, but on its magnitude k = oj/c only. The two terms 
/o and fi determine, respectively, the photon number 
density n — pfo and flux density j = ^cpii, where 
p{uj) — 4:Ttlu'^{2ttc)~^ is the density of states. Integra- 
tion of Eq. (I) gives two relations between n and i, 



^•j = ^C'(P/cq-")+PA, 



(7) 
(8) 



^c^r and mean free 



where the diffusion constant D 
path I — CT are determined by the transport scattering 
rate r ^ ^ — J^k' ""^kk' ( 1 — k • k' ) ■ The absorption length 
is defined by D^^^ — — . (An amplifying medium 
has an imaginary and a negative /cq.) In Eq. (|^) we 
have neglected terms of order (l/^a)'^, which are assumed 
to be < 1. 

Both Eqs. (^ and (^) contain a fluctuating source 
term. These two terms Cq and Ci have zero mean and 
correlators that follow from Eq. m), 
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(9b) 



where we have abbreviated A' — 6(lo — U!')6{t — t')S{r — 
r'). The correlator (9b) differs from the electronic case 
p3|-p^ by the factor 1 + fo instead of 1 — /o. This is the 
expected difference between boson and fermion statis- 
tics. The correlators (^^ and (^) have no electronic 
counterpart. They describe the statistics of the vacuum 
fluctuations. 

To demonstrate how the kinetic theory presented 
above works in a specific situation we consider the prop- 
agation through an absorbing or amplifying disordered 
waveguide (length L). The incident radiation is isotropic. 
All transmitted radiation is absorbed by a photodetector 
(see Fig. 1). Because of the one-dimensionality of the 
geometry we need to consider only the a;— dependence of 
j and n (we assume a unit cross-sectional area). The 
transmitted photon flux I = duj j{uj,L,t) fluctuates 
around its time-averaged value, I{t) — I + 61 (t). The 
(zero-frequency) noise power P — dt SI{t)5I{0) is 
the correlator of the fluctuating flux. We will compute 
P by solving the differential equations (^ and (^) with 
boundary conditions 71(0;, 0,t) — nin(a;,i), n{uj,L,t) — 0, 
dictated by the incident radiation at one end of the 
waveguide and by the absorbing photodetector at the 
other end. 




x = 



x = L 



FIG. 1. Isotropic radiation (solid arrows) is incident on 
a waveguide containing an absorbing or amplifying random 
medium. The transmitted radiation (dashed arrows) is ab- 
sorbed by a photodetector. 

Combining Eqs. (|^) and (||) we find equations for the 
mean and the fluctuations of the photon number. 
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(PSn Sn p dCi pCo 
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(10) 
(11) 



The homogeneous differential equation has Green func- 
tion 

r.( c sinh(a;</Ca)sinh(s-a;>/^a) 

U{X,X } — —t;a r-r , U^J 

smh s 
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where we have defined s = L/^^a and (a^>) is the 
smallest (largest) of x and x' . (In the amplifying sys- 
tem is imaginary so the hyperbolic functions become 
trigonometric functions.) The inhomogeneous equations 
( p^ ) and (|l^) have the solution 



n(cj, x) ~ . ( sinh s — sinh (x/£^a) — sinh (s — x/£,a 
smh s V 



+ ?^in(w) 



Pfc. 

iinh 

sinh(s - x/^a) 

sinh s 
L 



(13) 



Sn{uj,x,t) = p y"dx' G(a;,a:') - 



sinh (s - x/^g) 
sinh s 



(14) 



The flux density at the photodetector follows from Eq. 
(0) at X = L, 



Dm, 



sinh s ' 



(15) 



DSn\a Dp 



sinh s sinh s 



j dx [sin\i{x/£,a)-^ 



C0sh(x/^a)-^ 



(16) 



(Notice that the extra term cx Ci in Eq. (|^) is cancelled 
by the delta function in d'^G/dxdx' .) 

The time-averaged flux / = Jin -f /th is the sum 
of two contributions, the transmitted incident flux 
Jin = duj I?7iin/(Ca sinh s), and the thermal flux 
/th = ^^(Dp/cq/^a) tanh (s/2). The transmitted 
incident flux per frequency interval is a fraction T = 
4D/(c^Q sinhs) of the incident flux density jo = \cfiin- 
A fraction i? = 1 — 4£'/(c^a tanh s) of the incident flux 
is reflected. The thermal flux per frequency interval is a 
fraction 1-T - R ^ (4i:)/c^a) tanh (s/2) of the black- 
body flux density jo — jcpfcq- This is Kirchhoff's law of 
thermal radiation. 

The noise power P follows from the auto-correlators 
of £-0 and Ci [given by Eq. (^), with /o = n/p from 
Eq. (p^j. The auto-correlator of Sn-m and the cross- 
correlator of Co and Ci contribute only to order {l/^a)"^ 
and can therefore be neglected. The noise power P = 
Pin + Pth + Pex is fouud to cousist of three terms, given 
by 



-Pin = lin + du) 



Dn?^ 2s cosh (2s) + sinh (2s) - 4s 



sinh^ s 



Pth = Ith + dw 



Dpf^^ sinh^ (s/2) 
sinh'' s 



(17a) 
8s + 4s coshs 



7 sinh s — 4 sinh (2s) + sinh (3s) 



(17b) 



Per. 



, DfcqUin sinh (s/2) / 

did — 7 1 —OS — 4s cosh s 

2?a sinh s V 



+ 4 sinh s 4- 3 sinh (2s) 



(17c) 



The two terms Pin and Pth describe separately the fluctu- 
ations in the transmitted incident flux and in the thermal 
flux. Both terms are greater than the Poisson noise (the 
mean photon flux /th, /in) as a consequence of photon 
bunching. The third term Pex is the excess noise which in 
a quantum optical formulation originates from the beat- 
ing of the incident radiation with vacuum fluctuations in 
the medium [p^ . Here we find this excess noise from the 
semiclassical radiative transfer theory. The expressions 
for Pth and Pcx in Eq. (|l^) are the same as those that 
follow from the fully quantum optical treatment |pT| , |l7| . 
This is a crucial test of the validity of the semiclassical 
theory. The expression for Pin agrees with the quantum 
optical theory for the case that the incident radiation 
originates from a thermal source ||l^. The case of co- 
herent incident radiation is beyond the reach of radiative 
transfer theory. 

We envisage a variety of applications for the 
Boltzmann-Langevin equation for photons obtained in 
this paper. Although we have concentrated here on the 
waveguide geometry, in order to be able to compare with 
results in the literature, the calculation of the noise power 
in the diffusion approximation can be readily generalized 
to arbitrary geometry. As an example, we give the noise 
power of the thermal radiation emitted by a sphere (per 
unit surface area). 



Pth = /th + / duj / dz f cosh z 

J ta sinh s J 



sinh^ z 



(18) 



where s = P/^^ is the ratio of the radius R of the sphere 
and the absorption length ^q. The mean thermal flux is 
given by /th = duj /?p/oqC^(cotanhs - 1/s). The 
result for /th could have been obtained from the conven- 
tional radiative transfer theory using Kirchhoff's law, but 
the result for Pth could not. 

A dimensionless measure of the magnitude of the 
photon flux fluctuations is the Mandel parameter p9| , 
Q = [P — I)/ 1. In a photocoimt experiment, count- 
ing n photons in a time t with unit quantum efficiency, 
the Mandel parameter is obtained from the mean pho- 
tocount n and the variance var n in the long-time limit: 
Q = limt^oo(var n — n)/n. We assume a frequency- 
resolved measurement, so that the integrals over fre- 
quency in Eqs. (jl^) and (^8|) can be omitted. The Mandel 
parameter for thermal radiation from a waveguide and a 



3 



sphere is plotted in Fig. 2, as a function of s (s = L/^a 
for the waveguide and s = R/£,a for the sphere). Both 
the small and large-s behaviour of Q is geometry inde- 
pendent: Q = -^s^feq for s <C 1 and Q = ^feq for s ^ 1. 
The Bose-Einstein function feq{oj,T) is to be evaluated 
at the detection frequency ui and temperature T of the 
medium. The plot in Fig. 2 is for fcq = 10~^, typical for 
optical frequencies at 3000 K. 



extension to fluctuations presented here could be useful 
in that context as well. 
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FIG. 2. Mandel parameter Q ^ {P - I)/I for the thermal 
radiation from an absorbing medium and for the amplified 
spontaneous emission from a medium with a complete popu- 
lation inversion. The solid curves are for the sphere geometry 
[Eq. ( l^j , the dashed curves are for the waveguide geometry 



[Eq. (17t)]. The parameter s is the ratio of the radius of the 
sphere or of the length of the waveguide to the absorption or 
amplification length. The laser threshold in the amplifying 
case is at s = TT. To show both cases in one figure, the Q 
for the absorbing medium has been rescaled by a factor 10^ 



(corresponding to 



10" 



Much larger Mandel parameters can be obtained in 
amplifying systems, such as a random laser. Since com- 
plete population inversion corresponds to T — 0~, one 
has /eq = — 1 in that case 12 1. Eqs. ( |T7| ) and (|l8|) apply 



to amplified spontaneous emission below the laser thresh- 
old if one uses an imaginary ^a- The absolute value j^^l is 
the amplification length, and we denote s = L/\^a\ for the 
waveguide geometry and s — R/\^a\ for the sphere. The 
laser threshold occurs at s = tt in both geometries. We 
have included in Fig. 2 the Mandel parameter for these 
two amplifying systems for the case of complete popula- 
tion inversion. Again the result is geometry independent 
for small s, Q ~ ^s^|/cq| for s <^1. At the laser thresh- 
old (s = tt) the Mandel parameter diverges in the the- 
ory considered here. An important extension for future 
work is to include the non-linearities that become of cru- 
cial importance above the laser threshold. The simplic- 
ity of the radiative transfer theory developed here makes 
it a promising tool for the exploration of the non-linear 
regime in a random laser. 

Since radiative transfer theory was originally devel- 
oped for applications in astrophysics, we imagine that the 
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